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ABSTRACT 

We present a simple phenomenological model of feedback in early-type galaxies that tracks the evolution of the 
interstellar medium gas mass, metallicity, and temperature. Modeling the star formation rate as a Schmidt law with 
a temperature-dependent efficiency, we find that intermittent episodes of star formation are common in moderate- 
size ellipticals. Our model is applicable in the case in which the thermalization time from SN is sufficiently long 
that spatial variations are relatively unimportant, an appropriate assumption for the empirical parameters adopted 
here, but one that can only be demonstrated conclusively though more detailed numerical studies. The departure 
from a standard scenario of passive evolution implies significantly younger luminosity-weighted ages for the 
stellar populations of low-mass galaxies at moderate redshifts, even though the more physically meaningful mass- 
weighted ages are changed only slightly. Secondary bursts of star formation also lead to a natural explanation 
of the large scatter in the NUV-optical relation observed in clusters at moderate redshift and account for the 
population of E-i-A galaxies that display a spheroidal morphology. As the late-time formation of stars in our model 
is due to the gradual cooling of the interstellar medium, which is heated to temperatures ^ 1 keV by the initial burst 
of supernovae, our conclusions do not rely on any environmental effects or external mechanisms. Furthermore, a 
simple estimate of the X-ray emission from this supernova heated gas leads to an Lx vs Lb correlation that is in 
good agreement with observed values. Thus feedback processes may be essential to understanding the observed 
properties of early-type galaxies from the optical to the X-ray. 

Subject headings: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: stellar content 



L INTRODUCTION 

Of all the fields of astronomy, the study of galaxy formation 
may be the one in which the rift between theory and observa- 
tion is the most difficult to bridge. Theoretical astrophysicists, 
perhaps from a temperamental as well as a practical point of 
view, are best at predicting the most difficult component of the 
universe to measure: the overall distribution of invisible "dark 
matter" While the evolution of this component is well under- 
stood, numerical and analytical models must stretch their pre- 
dictive powers to the limit to superimpose on this history the 
messy heating, cooling, and enrichment processes that affect 
baryonic gas. 

At the same time, observational astrophysicists, perhaps 
from a romantic as well as practical point of view, are happi- 
est when looking at the most complicated thing they could ever 
measure: the stars in the night sky. A proverbial tail that wags 
the dog, the magnitudes and colors of the stellar populations 
in galaxies are dependent not only on the history of the bary- 
onic gas, but on the notoriously complicated process of star- 
formation. As many possible environmental effects can have an 
unknown impact on the number and distribution of stars formed 
(See, eg. Kroupa2001; Larson 1999; Scalo 1998) which is even 
more uncertain under primordial conditions (see, eg. Nakamura 
& Umemura 2001), predicting the optical properties of galaxies 
involved from "first principles" requires an enormous amount 
of extrapolation and simplifying assumptions. 

This great rift calls for astronomers to work to meet each 
other from each side of the divide. Thus simulations and semi- 
analytical models of galaxy formation are continuously striving 
to include all the important physical processes that affect form- 
ing galaxies and the surrounding intergalactic medium (e.g. 



Kauffmann et al. 1999; Baugh et al. 1998; Somerville & Pri- 
mack 1999). From the observational point of view, astronomers 
must constantly develop more careful phenomenological mod- 
els, which better express the features seen in simulations: mak- 
ing statements as to the ages of stars, their metallicities, and 
levels of dust extinction, rather than simply considering their 
infrared, optical, and ultraviolet colors and magnitudes. Mod- 
els of this sort include Tinsley (1980), Matteucci & Tornambe 
(1987), Tantalo et al. (1998), and Ferreras & Silk (2000b). 

From a theoretical point of view one of the most important 
issues to recently come to the fore is the role of supernova feed- 
back in galaxy formation. N-body simulations have shown that 
the baryonic components of galaxies that form in Cold Dark 
Matter (CDM) simulations systematically lose too much of 
their angular momentum to the dark matter halos in which they 
are contained (Navarro & Benz 1991). As discussed in White 
(1994), the resolution of this problem is widely believed to be 
the inclusion of feedback, although the proper modeling of this 
process remains a subject of intense debate. Thus since the first 
numerical model of thermal heating by Katz ( 1 992), approaches 
have ranged from implementation of kinetic boosts (Milos & 
Hernquist 1994; Gerritsen 1997), a multiphase model of star 
formation and feedback (Yepes 1997; Hultman & Pharasyn 
1999), a model of turbulent pressure support from feedback 
(Springel 2000), and a model in which energy persists in the 
interstellar medium for a period corresponding to the lifetime 
of stellar associations (Thacker & Couchman 2001). Finally, 
feedback processes are likely to play a more general role in im- 
pacting the intergalactic medium and conditions under which 
galaxy formation occurs (Mac Low & Ferrara 1999; Martel & 
Shapiro 2001; Scannapieco & Broadhurst 2001; Scannapieco, 
Thacker, & Davis 2001). 
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Yet despite the many exploratory feedback models being ex- 
amined theoretically, the role of feedback in empirical mod- 
eling has been more Umited. Ferreras & Silk (2000a; 2001) 
explored a simple model in which feedback was included phe- 
nomenologically using two free parameters: the star formation 
efficiency and the fraction of gas ejected in outflows. These 
values are dependent on the thermal and mechanical feedback 
from supemovae but they fail to capture their impact the state of 
the interstellar medium (ISM) itself, despite the close relation- 
ship between this state and star formation (McKee & Ostriker 
1977). Our primarily theoretical motivation in this paper, then, 
is to extend these models by tracking the temperature evolu- 
tion of the ISM, taking into account the thermal contributions 
of supemovae and metallicity-dependent gas cooling. 

The primarily observational motivation for this paper, on 
the other hand, is the measurement of a large scatter at the 
faint end of the near-ultraviolet (NUV) minus optical color- 
magnitude relation in cluster Abell 851 (z = 0.4) observed by 
Ferreras & Silk (2000a) using passbands F300W and F702W 
of the WFPC2 on board the Hubble Space Telescope. A sim- 
ple analysis of this scatter shows that it can not be accounted 
for by old Horizontal Branch stars. Moreover, significantly 
younger luminosity-weighted ages were found, although a de- 
generacy between the age of the young stellar component and 
its mass fraction prevented a detailed estimate of a more phys- 
ically meaningful mass-weighted age. This result shows that 
the star formation history of early-type galaxies is much more 
complicated than the standard picture in which a single stellar 
population formed at a redshift z ^ 3 and then evolved pas- 
sively. 

The structure of this work is as follows. In §2 we describe a 
simple thermal model of supernova feedback in elliptical galax- 
ies and its impact on the gas, stellar, and metallicity evolution 
of these objects. In §3 we describe the general features of our 
model, and the physics on which they depend. Our evolution- 
ary tracks are combined with population-synthesis models in 
§4 and used to study the color-magnitude relationship in ellip- 
ticals in both the optical and near-ultraviolet. In §5 we exam- 
ine the compatibility of our model with X-ray observations of 
early-type galaxies. In §6 we discuss the implications of our 
modeling for future observational studies of feedback in ellip- 
ticals as well as how our model compares to other theoretical 
approaches, and conclusions are given in §7. 

2. THERMAL MODEL OF FEEDBACK IN ELLIPTICAL GALAXIES 
2. 1 . Basic Assumptions 

In order to reconstruct the star-formation and chemical en- 
richment of early-type galaxies we follow an approach similar 
to that adopted in Ferreras & Silk (2000b) which in turn is an 
adaptation of the formalism of Tinsley (1980) in which the evo- 
lution within a galaxy is reduced to a few general parameters. 
Our model consists of only gas and stars. For each component 
we trace the net metallicity, counting all elements heavier than 
hefium in the same way. The mass in gas - Mg{t) - and in stars 
- Ms(t) - in each model are normalized to the initial total gas 
mass Mgo such that /ig(f) = Mg(t)/Mga and /^((f) = Ms(t)/Mg(). 

The gas component is fueled by infall of pre-enriched gas at 
a rate assumed to be a Gaussian function: 
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where t/ is the infall timescale and f/o is the epoch of maxi- 
mum infall, parameterized in terms of a formation redshift zf 
such that tizp) = tfo- We additionally assume a Salpeter (1955) 
Initial Mass Function (IMF) with cutoffs at 0. 1 and IOOM0, and 
an initial metallicity of the infall gas of Zo = Zq/10. This value 
is motivated by Renzini (1999) who used an estimated Zq/3 
metallicity of the z = Universe and the fact that ~ 30% of aU 
stars having formed at z ^ 3, to imply that the average metal- 
licity of the high-redshift Universe is~ l/3x 1/3^ 1/10 so- 
lar. Similar levels of pre-enrichment are consistent with the ob- 
served lack of G-dwarf stars with metallicities below ~ O.IZ© 
in the Milky Way disk, and the paucity of such stars in other 
massive early and late-type galaxies (Worthey, Dorman, & 
Jones 1996; Thomas, Greggio, & Bender 1999). 

2.2. Gas and Stellar Evolution 

We consider an infall model in which the star formation ef- 
ficiency is determined by the temperatures of the gas. In this 
case the limiting factor in star formation is not the conversion 
of cool gas into stars, but rather the cooling of the gas itself into 
molecular clouds. We write the star formation rate as a Schmidt 
law with respect to the normalized gas mass: 



(2) 



where Ceff(7') is the temperature-dependent star formation effi- 
ciency and 1 <J\f<2. A comparison of the Schmidt law with 
the observed star formation rate in local galaxies gives a value 
J\f 1.5 (Kennicutt 1998), adopted throughout this paper We 
assume a Fermi-step function for the dependence of the star 
formation efficiency on the temperature: 
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where T4 = T /lO'^K and a is a parameter that controls the steep- 
ness of the step. This function assumes that the star formation 
efficiency decreases abruptly for an average temperature of the 
ISM above T4 ^ 1, chosen to correspond to the temperature at 
which hydrogen is ionized. Finally, we assume fixed values 
for the normalization (C^g = 100), and the steepness of the step 
(a = 3). 

With these definitions the equations for the evolution of 
gas and stars can be written as a sum of gas accretion, star- 
formation, and ejection of material from stars: 



^ =/(0-Ceff(r4)Mf + (1-Bout)£(0, 
^=Ceff(T4)Mf-£(0, 



(4) 
(5) 



where the integral E{t) is the mass of gas ejected at a time t 
from stars at the end of their lifetimes: 



E{t). 
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dm(j)im){m - w„)V'(? - - Ts„) , 



(6) 



where (j)(m) is the initial mass function, Wm is the mass of a 
stellar remnant with an initial mass m, and nit is the mass corre- 
sponding to a stellar lifetime t, each of which can be fixed using 
the values as in Ferreras & Silk (2000b). A delay Tsn is included 
to account for the fact that it takes some time for the energy 
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from supernovae to spread throughout the galaxy. This should 
be roughly equal the time it takes for a gas bubble heated by the 
surrounding supernovae ejecta to percolate throughout the ISM 
of the galaxy, which we estimate to be t^^ ^100 Myr Finally, 
Bout is the fraction of gas ejected in outflows, as discussed in 
Ferreras & Silk (2000b). This parameter is a decreasing func- 
tion of galaxy mass that lies in the range < Bout < 1 and de- 
termines the average metallicity of the stellar populations. This 
inefficiency of metal ejection in large ellipticals is thought to 
be the best explanation of the correlation between mass and 
metallicity (Arimoto & Yoshii 1987) which gives rise to the 
color-magnitude relation. 

2.3. Thermal Feedback 

In order to model the temperature evolution of the gas we 
consider the cooling and heating sources that contribute to its 
total energy content. Defining 8 as the total internal energy, we 
find: 



^ = -N,n, K(T,Z) + ^ kTf+e lO^'erg 
dt Mq Mq 



(7) 



where Ng and He are the total number and number density of 
electrons, and Tf is the temperature of the infalling gas, and e 
is the fraction of the SN energy which goes toward heating the 
gas. The three terms on the right-hand side of this equation 
correspond to radiative cooling in the gas, thermal input from 
infalling gas, and energy input from SNe, respectively. 

The first of these is proportional to the electron density of the 
gas squared times the radiative cooling function, which is de- 
pendent on both the temperature and overall metallicity of the 
gas. In our simple model, we assume a polytropic equation so 
that rie oc r'/*'''"'\ The value of the polytropic index 7 depends 
on the state of the ISM as shown below. The cooling function, 
on the other hand, can be calculated directly according to the 
tabulated model in Sutherland & Dopita (1993), which assumes 
solar abundance ratios. 

In the second term, we take Tf to be the virial temperature 
of the halo; and in the last term, R^n is the rate of supernovae 
in the galaxy and e is the fraction of the resulting mechanical 
energy that is deposited into the gas. In calculating /?sn we con- 
sider contributions from both Type II and Type la SNe, such 
that Rsait) = Rii{t) + Ria{t)- The first of these contributions can 
be written as 



Rn{t) = I dm(j)(m)Ces[T4it-TmTsn)]fi-^it-T„r-Tsn), (8) 



where the integral gives the contribution from massive stars that 
undergo core collapse. 

In the Type la case, we follow the standard prescription of 
Greggio & Renzini (1983), recently explored in more detail by 
Matteucci & Recchi (2001). In this model the rate is written 
with respect to the distribution function of binary systems that 
can harbor a white dwarf: 

/> 1 6Mq /> ^iiicix 

Ria(t) = A dm(t)(m) dn'f(p')'llj(t-T„,-Tsn), (9) 

J mi J fi(m) 

where m, = max(3MQ , m(f )) and fj, is the ratio of the mass 
of the secondary star and the total, which ranges in the inte- 
gral from /i(m) = max(l -8/m,0.8/m) to = 0.5. Finally, 
/(/z) = 2'^''(1 + P)ijI^ is the distribution function of binaries 



which is taken from Greggio & Renzini (1983), with (3 = 2. 
The normalization in this case is taken to be ^ = 0.05, cali- 
brated from the ratio of Type la to Type II supernovae that best 
fits solar abundances N{Ia)/N{II) = Q. 12 (Nomoto, Iwamoto, & 
Kishimoto 1997). 

In order to compute the evolution of the temperature, we as- 
sume an ideal gas {p = nkT), with an equation of state given by 
p cx n'', which implies T cx n'''^ Writing 



d£ = -NkdT-pdV, 



we find 



d£ , /3 
— = dTi - + - 

kN \2 7 
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(10) 
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During the formation history of the galaxy we assume a sin- 
gle equation of state, namely, that of an ideal gas at constant 
pressure. Hence we fix 7 = throughout. This is motivated by 
cooling timescales which are much longer than the dynamical 
time. Rewriting the feedback efficiency in terms of the fraction 
of gas ejected in outflows, and solving for the temperature in 
terms of T4, the evolution of the gas temperature becomes 



2^?^ = (l-Bout)mQ.7Ts,„(Rn+Rj„)-[ip-(l-Bout)E]T4- 

-^,^+m-T4), (12) 



where mo 7 is the mean mass of a gas particle in units of 0.7 
times the proton mass (note that this gas is ionized), TcooI = 
kT/nA is the cooling timescale, and Ts„ ^ 3 x 10^ (in units 
lO'* K) is again the energy input from supernovae, now rescaled 
in terms of a supernova "temperature," which corresponds to a 
thermal efficiency of e = 0.7(1 -Bout)- 

Note that our model ignores the thermal contribution from 
the ejecta of intermediate and low-mass stars. This approxima- 
tion is based on the assumption that the temperatures of plane- 
tary nebulae are much less than the typical temperature of gas 
in the galaxy. 

2.4. Chemical Enrichment 

The final element of our model is the evolution of metals in 
the gas. In this case the relevant equation is 

=Z^/-ZjCeff(r)Aif +(l-1.2Bout)(£z+3'/«-R/„), (13) 

where Zf = 0. \Zq is the metallicity of the pre-enriched infalling 
gas. Note that each of these terms parallels those in eq. (Q). In 
eq. ([l3|) Ez is the amount of metals ejected from stars: 



Ez{t) 



dm4>(m)Tl>(t - Tm - Tsn) X 



X \(m-w„^ 



- mp„)Zg(t - T,„ - Tsn) + mp„] , 



(14) 



where he fraction of a star of mass m transformed into metals 
is given by p„,. For Intermediate Mass Stars (IMS), with M < 
8M0 we use a power law fit to the yields from Renzini & Voli 
(1981), Marigo, Bressan & Chiosi (1996) and van den Hoek 
& Groenewegen (1997). For larger stars, which contribute to 
metal enrichment via Type-II supernovae, the yields are taken 
from Woosley & Weaver (1995) and Thielemann, Nomoto & 
Hashimoto (1996). Note that the differences in the yields be- 
tween these two groups are due to the different physical inputs 
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assumed; mainly the criterion for convection, the nuclear reac- 
tion rates, and the determination of the upper mass cut. Finally, 
the (1-1 .2Bout) factor in eq. (|l3|) accounts for the fact that out- 
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Fig. 1. — Star formation history of a massive elliptical galaxy, which cor- 
responds to To = 300; no = 2; ry = 0.5 Gyr; zf = 2; Bom = 0. The photometry 
at zero redshift U -V = 1.55; V — K = 3.25 corresponds to a bright elliptical 
(Mv = —22; logo- = 2.4) in a local cluster. 



flows driven by supernova-induced winds are likely to be en- 
riched with metals, as described by Vader (1986). Finally, 
Type la supernovae must be added to the chemical enrichment 
process. Using the rate of supernovae - Ria{t) - as given by eq. 
(H), we include the contribution to the metallicity of the gas, as- 
suming the yield from each type la supernova to be yja = 0.6Mq 
(mostly iron; Thielemann, Nomoto & Yokoi 1986). There is no 
such term in the equation tracing the evolution of gas, eq. (^, 
however, since this type of supernova does not contribute to the 
gas reservoir. The standard progenitor of a type la supernova is 
a binary system with a C-O white dwarf. Given that the maxi- 
mum initial stellar mass that leads to a white dwarf is approxi- 
mately 8M0, this means type la supernova begin to contribute 
0.03 Gyr after stars are first formed. 

2.5. The Formation Epoch of Ellipticals 

In keeping with our desire to incorporate both observational 
and theoretical advances into our simple model, we do not ap- 
ply the usual observational benchmark of a fixed formation 
epoch, but instead allow for a range of formation redshifts cor- 
responding to the theoretically favored model of hierarchical 
structure formation. In this case we allow the formation redshift 
to be a function of the virial temperature, and choose zpiTwii) 
such that the probability of forming a galaxy at zf is the same 
for all values of Tvi, . Throughout this paper, the formation red- 
shift corresponds to the epoch of maximum gas infall. 

We restrict our attention to the currently favored ACDM 
model of structure formation, which is described in detail in 
Eisenstein & Hu (1999). In this model, objects that are equally 
likely to form have the same value of (t{R)D{zf), where (j{R) is 
the amplitude of mass fluctuations inside a sphere of radius R 



extrapolated linearly to z = 0, given by their eq. (34) and (35), 
and D is a "growth factor" which tracks the linear growth of 
perturbations as a function of redshift, given by their eq. (10). 
Based mainly on the most recent measurements of the Cosmic 
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Fig. 2. — Star formation history of a faint elliptical galaxy, which corre- 
sponds to To = 45; no = 0.8; Tf = 0.2 Gyr; zf = 3.9; Bout = 0.5. The photometry 
at zero redshift U - V = 1.2, V — K = 2.9 corresponds to a faint My = -18 elUp- 
tical (i.e. log cr ~ 1.8) in a local cluster. 



Microwave Background (Balbi et al. 2000; Netterfield et al. 
2001; Pryke et al. 2002), we consider a cosmological model 
in which n„, = 0.35, = 0.65, n,, = 0.06, erg = 0.87, T = 0. 18, 
n = 1, and /1 = 0.65, where il,„, Ha, and flh are the total matter, 
vacuum, and baryonic densities in units of the critical density, 
(78 is a(r) at the 8/i~' Mpc scale, h is the Hubble constant in 
units of 100 km s"' Mpc"', and n is the "tilt" of the primordial 
power spectrum, where n = 1 corresponds to a scale-invariant 
spectrum. The age of the universe in this cosmology is ^ 14 
Gyr. 

Using this model, we select a range of objects and determine 
Zf and T4 parametric ally as a function of total mass, which is not 
itself observable. To calculate zf, we compute the radius of the 
perturbation, R, as 0.95r2m '^^mJ^^ /z"' Mpc, where M12 is the 
mass in units of 1O'^/;~'M0, and then invert a(R)D(zF) = const. 
Similarly, we calculate from a standard model of collapse 
and virialization, which has been shown to be in good agree- 
ment with more detailed numerical simulations (See e.g.. Eke, 
Cole & Frenk 1996). In this case 



= 90KMif (1 +zf) [a„ + r!A(l +ZFr'] 



3nI/3 



(15) 



Finally, we select a(R)D(zF) such that the largest objects we 
consider — for which T4 = 300 — form at a redshift of zf = 2, 
leading to a range of formation redshifts from 2 to 4.4 (at 
T4 = 30). Note that in this case ellipticals are formed from the 
so-called "rare peaks" at redshifts which the mean amplitude 
of mass fluctuations, a(R)D(zF), is 2.5 times smaller than the 
critical value for collapse and virialization. This is consistent 
with hierarchical models which predict ellipticals to form at the 
peaks of the density distribution due to their clustering proper- 
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Table 1 

Mass-independent model parameters 



Parameter 


Symbol 


Value 


Schmidt-law index 


Af 


1.5 


Pre-enrichment 


Zo 


Zq/10 


SN Temperature 


TsmA 


3 X 10^ 


SN Delay 




100 Myr 


SF Efficiency 


^eff 


100 


Efficiency Steepness 


a 


3 



ties and the fact that such perturbations are likely to undergo a 
"major merger" (see eg., Kauffmann & Chariot 1998). 



3. INTERMITTENT STAR-FORMATION IN ELLIPTICAL GALAXIES 

3.1. Global Parameters and Features of Feedback Models 

In this section, we describe the main features of our model 
and the physics on which they depend. In Table 1 we list the 
model parameters that are independent of mass. The first two 
of these Af = 1.5 and Zq = Zq/10 are set to canonical values 
consistent with a variety of observations. The supernovae pa- 
rameters rsn.4 and Tsn, on the other hand, are based on simple es- 
timates of the thermal energy released by each SNe and the time 
delay for this energy to be transferred to the interstellar medium 
of the galaxy. Finally the thermal SF efficiency C^^f and steep- 
ness a are arbitrary parameters in our model, whose values we 
have fixed after an extensive search of parameter space. 

Having established these values, we show in Figure |l] the 
model prediction for the evolution of a massive elliptical galaxy 
(no = 2; Bout = 0; To = 300; zf = 2). This object undergoes a num- 
ber of bursts of star formation at roughly constant time intervals 
from its formation redshift of z/r = 2 until the temperature of the 
ISM rises above ^ 1 keV. Notice that the feedback cycles gener- 
ated by our model are quite punctuated, with the gas remaining 
at high temperatures for the majority of their evolution, inter- 
rupted by brief intervals of "catastrophic" cooling that result in 
new bursts of star formation. Each burst is followed by a new 
generation of supernovae, which quickly heats the ISM back up 
to the high temperature state, where it remains quiescent until 
the next epoch of cooling. 

Each phase of this feedback cycle is dependent on the physi- 
cal parameters of our models. In most cases, this dependence is 
quite straightforward and can be understood by simple scaling 
arguments, despite the fact that a more complicated numeri- 
cal approach is required to study the evolution in detail. Be- 
fore turning our attention to the observable consequences of 
our modeling, we first review these dependencies, comparing 
the properties of the large galaxy shown in Figure |l] with an 
example of a much smaller object with parameters (no = 0.8; 
Bout = 0.5; To = 45; zf = 3.9) in Figure | 

Both of these galaxies start at their virial temperatures, and 
remain there for a significant period of time, before undergoing 
a first burst of star-formation. At high temperatures ^ 10^ K, 
the cooling function A roughly scales as l/T, thus the cool- 
ing time of the gas TcooI = kT /nA is approximately cx T"-', as 
n (X l/T during the cooling phase. This strong dependence 



on temperature means that the main barrier to star-formation 
in these objects is the cooling time at the highest temperature 
in the feedback cycle. Once the galaxy cools to temperatures 
T ~ lO^K, the sharp decline of the cooling function halts the 
thermal run-away, while at the same time star-formation be- 
comes much more efficient. Note that the strong scaling of the 
cooling function with temperature means that this lower value 
is nearly constant over the full range of galaxy masses. It is 
during this stage that almost all stars are formed, resulting in 
a starburst with total star formation rates of order 1000 Mq/yt 
in massive galaxies and 10 Mo/yr in low-mass systems. The 
galaxy remains in a low-temperature state until the first genera- 
tion of supernovae explode and transfer their thermal energy 
to the gas. Thus each starburst results in a total fraction of 
stars created, which is roughly Ceff x (rs„ + T„,=iooM0) x figitss), 
where /i^(fsB) is gas mass at the start of the starburst. Simi- 
larly, the post starburst temperature is simply proportional to 
the number of stars formed times the energy transferred to the 

gas: £ CX (1 -BoudT^n X Ceff X (rsn+r„=i00Mo) X ^J.g(tsB). 

From these scaUngs, it is easy to understand why smaller 
galaxies oscillate more frequently than those with larger 
masses. As the energy per supernova, star formation efficiency, 
and Tsn are fixed for all models, the supernova energy generated 
in each model varies only by a factor (1 -Bout) X l^gitss). But for 
each starburst, the ratio of the initial and final gas mass is fixed 
and the post starburst temperature scales simply as cx (1 - Bout)- 
Thus the period between starbursts can be estimated as approx- 
imately cx r-' cx (1 -Bout)^. 

Comparing Figures 1 and 2, which give the star formation 
histories of two early-type galaxies with very different masses, 
we see that the ratio of gas outflow fractions is 2, implying that 
the final bursts of star formation will take place approximately 
8 times later in the high-mass galaxy. Thus the last early burst 
in the large galaxy takes place after a delay of approximately 1 
Gyr, while the last early burst in the T4 = 45 model is partially 
blended with the peak before it, with a relative separation of ap- 
proximately 1/8 Gyr Similarly, after the early burst stage, the 
smaller galaxy reaches a temperature of 1.5 x lO^K, and cools 
to form a secondary late burst of star formation after approxi- 
mately 7 Gyrs. On the other hand, the large galaxy reaches a 
temperature of « 3 x lO^K after the initial burst phase, corre- 
sponding to a cooling time of ^ 50 Gyrs, and therefore no late 
bursts appear in this model. 

It is important to note that the relationship between mass and 
metallicity in our models arises not only from a differential Bout 
as in Ferreras and Silk (2000b), but from an interplay between 
the multiple bursts and the infall timescale. This aspect of our 
model is independent of the late time star formation properties 
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of these galaxies. Thus it is conceivable that the early time star- 
formation history of the galaxy may have been somewhat differ- 
ent (perhaps due to the merger of two disks for example) while 
still preserving the star formation cycles and scalings seen in 
our models at late times. The metallicity-mass relation is dis- 



o 




Fig. 3. — Average V-band luminosity-weighted age and metallicity (left) and 
its dilference witli a mass-weighted average (right) at zero redshift. The lines 
in the left panels give a least squares fit to the age-mass and metallicity-mass 
relations. The slopes of these Hnear fits are 0.02 and 0.12, respectively. 



cussed in more detail in §4. 

Finally, in Figures 1 and 2, the thick lines represent an es- 
timate as to the physical star formation rate in these objects. 
This is likely to be a smoothed version of the theoretical curve, 
as our model assumes a point-like galaxy, whereas a more re- 
alistic model including spatial information would have bursts 
of star formation happening at slightly different epochs in dif- 
ferent regions of the galaxy. This would result in a more con- 
tinuous global star formation rate during the first stage, while 
preserving the punctuated bursts at later times. 

This raises the question as to whether the interactions be- 
tween areas of the galaxy in which star formation is out of phase 
may be able to suppress the cycles seen in our one-zone mod- 
els. Here the important issue is whether SNe are able to heat 
a region of space and suppress star formation fast enough that 
spatial variations between regions become important. In our 
one-zone models, we assume a typical delay of 100 Myrs for 
SN to percolate through the galaxy, dumping a fraction Bout of 
their energy into the ISM. Thus we expect a one-zone descrip- 
tion to be applicable if the galaxy is relatively homogeneous 
on scales comparable to 100 Myrs times the sound speed of the 

1/2 

gas. A simple estimate of this speed gives w lOT^ km/s, our 



one-zone approach is equivalent to assuming that the galaxy is 

1/2 

relatively smooth on scales above tsnc, ~ IT^ kpc, a reason- 
ably large scale even at 74 = 1 . Note however that our value for 
the SNe delay time is somewhat arbitrary and chosen empiri- 
cally to reproduce the observed properties of elliptical galaxies. 
Thus while our model is a plausible description of a more com- 
plete treatment, a dedicated numerical study would be neces- 



sary to more firmly establish the relevant thermaUzation time in 
our models, and whether c^r is large enough that spatial varia- 
tions can be ignored. 



3.2. Mass Dependent Parameters and Average Properties 
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Fig. 4. — Same as figure g for a moderate redshift (z = 0.5) cluster. The 
slopes of the linear fits to the age-mass and metallicity-mass relations are 0.02 
and 0.20, respectively. 

Having examined the main features of our models, we now 
turn our attention to constructing a representative sample of el- 
liptical galaxies over a range of masses. In Table || we list the 
parameters of our model that depend on age and indicate their 
scaling with mass and formation redshift. 

The virial temperature Tvir simply scales as GM^I^p^l^, 
where p is the mean cosmological gas density at the time of 
formation. Similarly the initial number density of the gas at 
Kf'K, no oc pTyir as n oc l/T during the cooling phase. The in- 
fall timescale, on the other hand, scales as the radius of the halo 
oiM^^^p^^^ divided by the virial velocity cx T^l^. Finally, we as- 
sume a simple power law dependence with mass for the fraction 
of gas ejected in outflows, Bout- The most massive galaxies are 
expected to have negligible outflows (Bom ^ 0), while the low 
mass galaxies should have higher outflow fractions that will re- 
duce the average metallicity of the stellar populations, thereby 
blueing their colors. We use the color of faint ellipticals in local 
clusters to constrain Bout in these systems, predicting roughly 
Bout ~ 0.6 for Mio ^ 5. Thus the power law index over the 
mass range 5 < Mio < 300 is ~ -0.3. 

Usin g the mass-independent parameters described in Ta- 
ble 2.5 as well as the scaling of the mass-dependent parameters 
shown in Table ^, we can generate a sample of galaxies start- 
ing from a luminosity function. We estimate the relative num- 
ber densities of these objects using the NIR //-band luminosity 
function of the Coma cluster (De Propris et al. 1998), which 
provides a good estimate of the underlying mass function, and 
arrive at a final sample of galaxies studied in detail below. 

Figures || and ^ show mass and V-band luminosity-weighted 
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averaged ages and metallicities of the objects, defined as: 

rt(z) . rt(z) 

<tM>= / dttipit) I dti){t), (16) 

^0 ' JO 



<Zm> 



i{z) 



dtZip(t)^ [ dttpit) 
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-20 



(17) 
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relation is nearly flat both at zero and moderate (z = 0.5) red- 
shift (f oc M"-"^), although the latter presents a significantly 
larger scatter. While there is little difference between mass and 
luminosity-weighted ages at zero redshift, Figure ^ shows sig- 
nificant differences at z = 0.5. The mass-metallicity correlation 
in this figure is clearly defined and has a similar power law in- 
dex: Z oc Af*'^^. However, there is an increased scatter in the 
luminosity-weighted ages and metallicities, especially at 
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Fig. 5. — Color-magnitude relation in optical and near-infrared passbands 
for a local cluster. The filled squai'es are our model predictions. The hollow 
circles are Coma cluster ellipticals observed by Bower et al. (1992). The best 
linear fit to the observations is given by the solid line, and the ±0.05 mag scatter 
found in these galaxies is plotted as a dashed line. 



/•t(z) . pt(z) 

<tL>= J dttT(t)-^il>(t) I J dtTity^iPit), (18) 

rt(z) rt(z) 

<Zl>= J dtZT{t)-^^{t) I jf t/fT(f)"V(0, (19) 

where T(f) is the V-band mass-to-light ratio as a function of 
stellar age, and t{z) is the age of the universe at the redshift at 
which it is observered. Comparing these results, it is clear that 
the ages of our model galaxies at moderate redshift are strongly 
dependent on which type of average is applied. While the mass- 
weighted values provide the most meaningful summary of the 
overall star-formation history of the object, both photometric 
and spectral observations only give estimates of the luminosity- 
weighted age and metallicity. Thus, late episodes of star forma- 
tion skew the integrated spectral energy distribution to younger 
ages, because of the very small mass-to-light ratios of high- 
mass main sequence stars. 

Figure |^ shows that at zero redshift, early-type galaxies can 
be described by a metallicity sequence roughly between 1.5Z0 
for the most massive galaxies, down to Q.5Zq for the faintest 
ones, in agreement with observations of ellipticals in local clus- 
ters (Kuntschner & Davies 1998; Kobayashi & Arimoto 1999). 
Our model gives a power law fit between mass and metallic- 
ity of Z ocM^ '^ at zero redshift. The slope of the age-mass 



Fig. 6. — Rest frame U-V (left) and NUV— optical (right) color-magnitude 
relations predicted at redshifts of 7, = 0.55 and z = 0.41. The filled squares are 
the model predictions, whereas the hollow circles show the data from cluster 
Abell 851 (z = 0.41) explored in the NUV by Ferreras & Silk (2000a) and 
C10016+16 (z = 0.55) observed by the MORPHS collaboration (Ellis et al. 
1997). The shaded areas in the panels on the right give the 4(7 detection limit 
from the shallower F300W images. The solid line in the panels on the right 
gives the prediction of a single-age metallicity sequence formed at Zf = 10 us- 
ing the color-magnitude of the Coma cluster as constraint. The solid and dashed 
hues in the bottom left panel is the fit and scatter, respectively, to the color mag- 
nitude relation in the elliptical galaxies of cluster C10016+16. 

the faint end. This offset can be as large as a few Gyr in age, es- 
pecially for galaxies with secondary bursts occurring within the 
previous 1 Gyr. This occurs since even small amounts of young 
stars can boost the luminosity mainly in the NUV, thereby re- 
ducing the luminosity-weighted age. 

4. OPTICAL AND NEAR-ULTRAVIOLET COLOR-MAGNITUDE 
RELATIONS IN ELLIPTICALS 

For a given choice of parameters, our model predicts a chem- 
ical enrichment track that determines the ages and metallicities 
of the stellar populations. These star formation histories can 
be convolved with models of simple stellar populations to get 
the spectrophotometric properties of the model galaxies. To do 
this, we use the latest population synthesis models of Bruzual & 
Chariot (in preparation) and assume a fixed baryon to total mass 
ratio of ilfe/il,„ ^ 1/6 in order to obtain absolute luminosities. 

Figure gives the optical and near-infrared (NIR) color- 
magnitude relation predicted for a local cluster The model 
galaxies are shown as filled squares and are compared with the 
precision photometry of Coma ellipticals observed by Bower, 
Lucey & Ellis (1992), plotted as hollow circles. The solid and 
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Table 2 

Mass-dependent model parameters 



Parameter 


Symbol 


Scaling 


Virial temperature 


7vir,4 


cxM^^'(l + Zf) 


Initial number density at r4 = 1 


no 


CX 7;ir,4(l+ZF)^ 


Infall timescale 




(x(l+ZFr3/2 


Outflow gas fraction 


Bout 


cx-0.31ogMio 



dashed lines are the best fits and scatter to the observations, 
respectively. The remarkable agreement shows that our Bout 
scaling reproduces well the expected mass-metallicity or mass- 
luminosity relationship (as shown in figure |]) that explains the 
color range in early-type galaxies. Notice that the V -K color 
departs towards the blue at the faint end more than the U -V 
color. This has to do with the higher age sensitivity of optical 
colors with respect to near-infrared passbands. The U and V- 
bands will be more affected by recent episodes of star formation 
than ^T-band, so that V-K appears bluer than U-V at the faint 
end at which late bursts of star formation occur. 

Figure || shows the rest-frame 2000A near-ultraviolet (NUV) 
and optical color-magnitude relation predicted for two clusters 
at moderate redshift: Abell 851 (z = 0.4; right) whose optical 
and NUV color-magnitude relation was recently analyzed by 
Ferreras & Silk (2000a), and C10016H-16 (z = 0.55; Jeff) whose 
detailed optical color-magnitude relation can be found in Ellis 
et al. (1997). The shaded area gives the 4(7 detection limit for 
the shallower F300W images. The solid and dashed lines show 
the fit and scatter to the observed data in cluster C10016H-16. 

These two color-magnitude relations paint two very different 
pictures of a similar cluster at moderate redshift. The optical 
colors shown in the panels on the left do not present a large 
scatter This is consistent with the observations of z < 1 clus- 
ters by Stanford, Eisenhardt, & Dickinson (1998), which they 
interpreted as an absence of star formation at redshifts below 
2-3. 

The near-ultraviolet color-magnitude relation showed in the 
panels on the right, however, point to a more involved SF his- 
tory. In this case the blue band maps a rest frame spectral win- 
dow around 2000A in which young main sequence A-type stars 
contribute significantly, making these colors much more sen- 
sitive to recent star formation. Thus the observed scatter with 
respect to the prediction for old stellar populations formed at 
redshifts z ~ 5 - 10 (thick line) is a telltale sign of recent star 
formation. This is particularly clear as the lookback time of this 
cluster allows us to safely ignore the contribution to the F300W 
flux from evolved, core helium-burning stars, whose contribu- 
tion would peak around 1500A. Hence the combination of NUV 
and optical colors of early-type galaxies at moderate redshifts is 
a valuable observable to infer the recent star formation histories 
of these systems. Yet detailed analyses of NUV-optical color- 
magnitude relations at moderate redshift are scarce (Buson et 
al. 2000; Ferreras & Silk 2000a) and more work is needed in 
this direction. 

Mass-to-light ratios are also useful stellar clocks although 
their observed estimates are highly model dependent. Fig- 
ure^ shows the model predictions against total mass in units of 
1O'*W0. Galaxies with late bursts depart from the "red enve- 



lope" since their young stars contribute significantly to the total 
V-band luminosity, even if they make up only a small fraction 
of the total stellar mass. Here the solid line is the corrected 
fit to the observed data from Mobasher et al. (1999; stars) and 
Pahre (1999; hollow squares) which gives M/Ly oc M^-^*. The 
correction involves a slightly nonlinear power law dependence 
between the stellar mass Ms and the total mass M oc M] ^ as 
described in Ferreras & Silk (2000b). 

5 . X-RAY EMISSION FROM EARLY-TYPE GALAXIES 

The model presented in this paper gives the thermal state of 
the ISM a central role in the star formation history of elliptical 
galaxies. Even though our model is a rather simple 1-zone sys- 
tem, we can use this information to predict the X-ray luminos- 
ity from the gas and its correlation with the optical luminosity, 
which traces the stellar component. 

The X-ray luminosity of the galaxy can be approximated 
as arising purely from Bremsstrahlung emission from free-free 
electron interactions. In this case 

Lx « 2 X lO~'^\,NeTl'^ [e'^'^ -e"T^^ ergscms"' , (20) 

where Zsmax and Emin are the maximum and minimum energies 
of the observed X-ray band. Due to line emission, this equation 
is only valid at temperatures much greater than r4 = 1, but this 
is a good assumption as the X-ray emission from cold objects is 
negligible. A larger source of error, in fact, is that we only in- 
clude X-ray emission from the gas and not from X-ray binaries, 
which may contribute significantly. 

Rewriting the number density in terms of temperature and n, 
using A^j. = 1 .5 X lO^^Mio in this cosmology, and restricting our 
attention to the 0.5 to 2.4 kev energy band, we find 

il 1 1/9 / -522 -3400 \ , 

Lx«3x lO'^'MioncmVjT/ fe ^4 -e ^4 jergss"'. (21) 

The filled squares in Figure ^ show the predicted Lx vs Lb cor- 
relation in early-type galaxies following this equation. The hol- 
low circles are taken from O'Sullivan et al. (2001) who com- 
piled a sample of 401 early-type galaxies. The solid line gives 
the best fit to the observed data, excluding bright cluster and 
group galaxies as well as AGNs. Our model is in good agree- 
ment, showing that the trend of the Lx '■ Lb correlation in ellip- 
tical galaxies can be mostly explained by X-ray emission from 
the hot interstellar medium. While the overall scatter in the data 
exceeds that in the models somewhat, this is not surprising, as 
the observational values are subject to additional variations not 
only due to the X-ray binaries as described above, but perhaps 
also due to environmental effects (Mathews & Brighenti 1998). 
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We can also use a simple scaling argument to estimate sizes 
from the galaxy mass (Mio) and number density (no « nT4). In 
this case 



- 5 kpc 



«()mo,7 



1/3 



(22) 



For our model galaxies the inferred sizes range from 8 to 
30 kpc, which is a reasonable estimate of the core X-ray sizes of 
ellipticals (Mathews & Brighenti 1998). High-resolution deep 
X-ray images of ellipticals over a large mass range would be 
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Fig. 7. — Correlation between the stellar V-band mass-to-light ratio as a 
function of total mass in units of IO'^Mq . The model data (filled squares) have 
been coirected for the offset between the total and the stellar mass content in 
elliptical galaxies (Fen'eras & Silk 2000b). The observations are from Pahre 
(1999; P99) and Mobasher et al. (1999; M99). 



needed to determine the scaling of galaxy X-ray size with opti- 
cal luminosity, which could quantify the contribution to the X- 
ray brightness from the hot interstellar gas in elUptical galaxies. 



6. TOWARD A MORE COMPLETE PICTURE OF FEEDBACK ESf 
EARLY-TYPE GALAXIES 

While our simple model of thermal feedback yields a promis- 
ing explanation for many of the observed properties of early- 
type galaxies, it naturally raises two questions. As secondary 
starbursts may also arise from more complicated gas accretion 
histories that the one consider here, one is left with the ques- 
tion of the what types of observations will be most important 
in distinguishing between our scenario and other possibilities. 
Secondly, from a theoretical point of view, it is clear that intense 
bursts of supernovae will deposit both energy and momentum 
over a much larger area than the ISM of the galaxy itself. This 
raises the question as to the uniqueness of our single zone ther- 
mal model, and whether similar feedback loops can also arise 
in the interplay between the galaxy and its environment. In this 
section we address each of these questions in turn. 

6.1. Observational Implications 



The elliptical SF histories described in this paper are much 
more complicated than those of more usual models which as- 
sume a short epoch of star formation resulting in stellar popu- 
lations with well-defined ages and metallicities. Although our 
predictions for massive galaxies are quite similar to such mod- 
els, our lower mass systems experience secondary episodes of 
star formation that affect their luminosities, particularly at mod- 
erate redshifts. This is true independent of the possibility of 
other sources of late star formation from environmental effects. 

Our conclusion clearly calls for more detailed analyses of el- 
liptical galaxies at moderate and high redshift. It also gives an 
additional mechanism for "hiding" high-redshift ellipticals 
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Fig. 8. — Predicted correlation between the X-ray and optical luminosi- 
ties using a simple argument involving Bremsstrahlung emission from the elec- 
trons in the hot gas component. The filled squares give our model predictions, 
whereas the hollow circles are from O'SuUivan et al. (2001). 



from view. In fact, this possibility is reminiscent of the the 
blue-nucleated field early-type galaxies observed by Menan- 
teau, Abraham & Ellis (2001) in the Hubble Deep Field, which 
may be examples of the late bursts predicted by our model. 
Similarly, the post-starburst E-i-A galaxies found in clusters are 
suggestive of our results. 

It is also interesting to note that the observed correlation be- 
tween abundance ratio enhancements and galaxy mass or lu- 
minosity (Kuntschner & Davies 1998; Trager et al. 2000) is 
indicative of extended star formation in low mass galaxies, 
which allows the byproducts from type la supernovae to be 
locked into subsequent generations of stars. Galaxies with cen- 
tral velocity dispersions around 50 km s"' , for example, tend 
to have [Mg/Fe] ratios close to Solar, while values around 
[Mg/Fe]^ -1-0.3 are observed in massive (ctq ^ 400 km s"') 
galaxies. Again, this hints at late bursts of star formation. 

While the short duration of our secondary bursts makes their 
direct detection a complicated endeavor, observations of NUV 
colors may provide a tractable approach to studying feedback. 
As young main-sequence A-type stars have strong Balmer ab- 
sorption lines in the NUV spectral window around 2000A, ob- 
servations that map this rest-frame wavelength can used to trace 
recent star formation. Thus while low numbers of cluster galax- 
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ies with strong emission lines are seen at moderate redshift 
(Dressier et al. 1999), up to 20% of such galaxies feature post- 
starburst NUV spectra (Poggianti et al. 1999). The large scatter 
at the faint end of the NUV minus optical color-magnitude rela- 
tion of Abell 85 1 (Ferreras & Silk 2000a) is a further indication 
of the sensitivity of this approach. As the thermal feedback 
mechanism presented in this paper is independent of the envi- 
ronment, comparisons of these colors in the field, groups, and 
cluster populations should provide a direct method of separat- 
ing the contribution of thermal feedback to star formation in 
spheroids from other processes. 

6.2. Alternative Models of Feedback 

The model of feedback in elliptical galaxies we have exam- 
ined in this paper focuses on ISM heating, modulating star for- 
mation by changing the condition of the gas within the galaxy 
itself. Yet the presence of a large number of supernovae fol- 
lowing a burst of star-formation is likely to also deposit sig- 
nificant amounts of energy and momentum into the surround- 
ing intergalactic medium. In our thermal model, infall onto the 
galaxy was fixed simply by the dynamical timescale of a viri- 
alized cloud as a function of mass and redshift. While this is a 
reasonable approximation at early times, the rate of infall from 
the intergalactic medium after the onset of star-formation may 
be slowed due to collisions with supernova ejecta. This raises 
the question as to whether multiple starbursts can also arise in 
models that focus on the interplay between the galaxy and its 
environment. 

To study the features of such models, we can replace our 
simple two-component model with a three-component model in 
which we track stars, the gas contained within the galaxy, and 
a second gas reservoir that represents the surrounding medium, 
which is condensing onto the galaxy. In order to isolate the 
features of such an approach from the thermal effects discussed 
above, we do not attempt to track the temperature of the gas, but 
rather consider a generic set of equations describing a broad 
class of models in which infall is modulated by feedback. In 
this case we find 

fig = flr/T-en,omE-Ce{ffJ,g + (l-f)(l-Boat)E, 
i^r = -fir/r+emomE + fil-BoudE, 

fis = Ceff/ig, (23) 



where emom is a parameter of order -y/ lO^'ergs/MQ/vvh- ~ 
10/vioo that parameterizes the suppression of infall by super- 
novae, / is the fraction of ejected material that is recycled into 
the infalling reservoir, and the infall is now parameterized as an 
exponential with a scale time t in order to simplify the analyti- 
cal arguments below. 

With these equations in hand we now approximate the be- 
havior of this system using a "post" instantaneous recycling ap- 
proximation in which we model E{t) considering the timescale 
(t2) corresponding to the lifetime of stars that play a role in 
stellar feedback. In this approximation we write the ejected 
fraction as 



£(?) ~ Ceff/sN 



l^git)-lj'g(t)T2 



(24) 



where /sn is the gas fraction returned by supernovae. Note that 
a good estimate of T2 is ~ 0.02 Gyr (roughly the lifetime of a 
massive star that undergoes core collapse, i.e. M ^ IOMq), but 
we are also free to take T2 = 0, recovering the more commonly 
studied instantaneous recycling approximation. 



Rewriting these equations in terms of the infall timescale and 
defining i = f /ti ; C = Ceff ti ; and x = T2 /ti , we find 

(l-CL]^-iif + (l+CM)^ + CNfi, = 0, (25) 
/ ds-^ ds 

where L = /sn e'x, M ee 1 + /sn e' + (l -Bout)/sN-^, = 1 - (1 - 
Bout)/sN, and e' = e^^m - (1 - /) (1 - ^out)- Note that for any rea- 
sonable choice of parameters e' > as emom ~ 10/vioo is much 
greater than (1 -/)(1 - Bout)- Thus L, M, and are > in all 
physical models. 

Eq. ( p5[ ) is the equation for a harmonic oscillator, y+'^y + 

WqJ = 0. If we define il = \J iJ^-^- /A-, then the criterion for 
oscillation is > 0, which can be rewritten as 

ACN{l-CL)-{l+CMf>Q. (26) 

In the limit in which C is very small or large this inequal- 
ity is never satisfied, and likewise at its local maximum C = 
(2N-M)/{M^+4NL). Thus we see that for very a general class 
of models in which infall is modulated by star formation, no 
oscillations are present in the post-instantaneous recycling ap- 
proximation. 

As an additional check, we have solved eq. ( |23| ) numeri- 
cally over a broad range of parameters, blanketing the relevant 
physical values. These results confirm the analytical arguments 
above, showing that oscillations never take place in models that 
modulate infall as long as e' = etnom-(l -/)(1 -Bout) is positive. 

7. CONCLUSIONS 

While stellar feedback in galaxy formation has been the sub- 
ject of intense theoretical investigation, these studies have had 
little impact on the direct interpretation of the observed prop- 
erties of elliptical galaxies. Yet there are several observational 
clues that point to its importance. Ferreras & Silk (2000b) for 
example, found that the color-magnitude relation of early-type 
cluster galaxies can be easily understood in the context of a 
variable ejection efficiency of supernova material, which scales 
inversely with galaxy mass. Similarly, the large scatter in the 
faint end of the near-ultraviolet minus optical color-magnitude 
relation in early-type cluster galaxies hints at episodic star- 
formation in these objects, a natural consequence of stellar 
feedback (FeiTeras & Silk 2000a). 

Motivated by these observational and theoretical pointers, we 
have explored in this work a simple model of thermal feed- 
back in elliptical galaxies. Our model consists of only gas and 
stars in a single zone and is a natural extension of previous 
investigations (Tinsely 1980; Ferreras & Silk 2000b; Ferreras 
& Silk 2001) which includes the thermal state of the interstel- 
lar medium. We account for heating by supernova feedback 
and line emission cooling, and their impact on the temperature 
and density of the ISM using a minimum of parameters, and 
relate these quantities to the overall star formation rate using 
a Schmidt law and a temperature-dependent efficiency. This 
model is applicable in the case in which the time for SNe to 
percolate through the ISM is sufficiently long that spatial vari- 
ations between regions are relatively unimportant. While this 
is the case for the choice of parameters adopted in this paper, 
more detailed modeling is necessary to study the interplay be- 
tween thermalization and substructure conclusively. 

Our model leads to two important and independent results. 
First the interplay between infall and variable SNe ejection ef- 
ficiency provides a natural explanation of the optical and NIR 
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color-magnitude relations in elliptical galaxies both locally and 
at moderate redshift (z <J 1). Our model generates a mass- 
metallicity relation Z oc m0.15-o.2^ which accounts for the ob- 
served color range over three magnitudes in luminosity. Fur- 
thermore, the rapid cooling times for temperatures close to 
10"* K result in short bursts of star formation, generating stellar 
populations that are very similar to the simple ones conmionly 
used to describe ellipticals. This implies no change in the slope 
or the scatter of the optical and NIR color-magnitude relation 
out to redshifts z ^ 1 as observed by many authors (e.g. Stan- 
ford et al. 1998; Van Dokkum et al. 1998; Van Dokkum et al. 
2000) 

Secondly, we find that secondary peaks of late star-formation 
are ubiquitous in smaller systems in which the fraction of the 
SNe ejecta that escapes from the ISM, Bout> is significant. In 
these cases the gas is heated to temperatures oc (1 -5out). lead- 
ing to secondary bursts of star formation with a delay that scales 
as oc (1 - -Bout)"' due to the scaling of the cooling times of 
high temperature gas. As such bursts occur within a Hubble 
time only in the smaller galaxies in which Bout is relatively 
large, this leads to a natural explanation of the large scatter 
in the NUV-optical relation as observed in clusters at moder- 
ate redshifts. This can also account for the observed population 
of post-starburst Eh-A galaxies that display a spheroidal mor- 
phology (Ferreras & Silk 2000c). While the current claim for 
these post-starburst systems is the quenching of star formation 
while falling through the hot intracluster medium (Poggianti et 
al. 1999) our model shows that late starbursts may arise in ellip- 
tical galaxies without resorting to environmental mechanisms. 
In fact, we have found that no such oscillations arise in a broad 
class of theoretical models which study the interplay SNe and 
the further accretion of gas. 

Given that the duration of these peaks of star formation is 
rather short (~ 100 Myr), the most promising observational ap- 
proach is to examine the NUV properties of spheroidal galaxies. 
A comprehensive study of the rest frame NUV-optical color- 
magnitude relation in ellipticals may be able to quantify both 
the number of that undergo late bursts and the mass fraction in 
young stars. Furthermore, as the thermal mechanism studied 
this paper is independent of environment, comparisons of such 
colors between field, group, and cluster populations can help 
to differentiate between this feedback and other processes. As 
the contribution in the NUV from evolved core helium-burning 
stars becomes significant with age, studies of galaxies at red- 
shifts z J> 0.3 may provide the cleanest samples for such com- 
parisons. 

At z ^ 1, late bursts can be significant even in high-mass 
galaxies, as shown in Figure 1. This may result in the failure 
of any search for high-redshift elUpticals that is based on sim- 



ple passively evolving models and does not invoke complicated 
scenarios of assembly (e.g. Zepf 1997). Our model also under- 
scores the bias intrinsic to observing luminosity weighted quan- 
tities. While the mass and luminosity weighted ages of higher 
mass ellipticals are quite similar, even the relatively small bursts 
of late star formation that arise in smaller objects cause large 
changes in their observed V band luminosity-weighted ages. 
This effect is even more severe in the NUV, highlighting the im- 
portance accounting for the overall star-formation history when 
interpreting observed stellar ages. 

One of the remarkable results of our model is the presence 
of an interstellar medium that has been heated by supemovae to 
temperatures around 1 keV or higher. This hot gas exists at very 
low densities, with only a fraction leaving the galaxy, much in 
the same way as the hot gas of the solar corona is kept gravitat- 
ing around the Sun. A full 3D simulation would be needed in 
order to explore the real conditions under which this "fountain 
effect" can occur. 

Finally our prediction of a hot corona of SNe heated gas has 
natural implications for the X-ray properties of elliptical galax- 
ies. By approximating the overall X-ray luminosity of each ob- 
ject as due to Bremsstrahlung emission from gas at the single 
temperature and density given by our model, we derive a Lx vs 
Lb correlation that is in good agreement with observed values. 
Nevertheless our simple one-zone model does not account for 
environmental effects and ignores what is likely to be a signif- 
icant contribution from unresolved X-ray binaries, resulting in 
a scatter that is smaller than observed. 

Galaxy formation is one of the richest and most complex 
processes in all of astrophysics, forcing observers and theo- 
rists to approach it from completely different viewpoints. And 
although both theoretical and observational progress has been 
clear and systematic, the rift between theory and observation 
in the field remains one of the most difficult to bridge. In this 
work, we have attempted to draw the key theoretical issue of 
feedback into the direct interpretation of the observational prop- 
erties of early-type galaxies. While our approach has been ex- 
plorational, the widespread agreement of our simple model with 
diverse observations suggests that thermal feedback processes 
are likely to be essential to fully understanding the optical, ul- 
traviolet, and X-ray properties of early-type galaxies. 
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